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ABSTRACT 

We have obtained high-resolution data of the z ~ 2 ring-like, clumpy star-forming galaxy (SFG) 
ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing- limited mode (in H- and 
J-band). Our data includes all of the main strong optical emission lines: [Oil], [OIII], Ha, Hf3, [Nil], 
and [SII] . We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's 
massive, star-forming clumps (<r ~ 85 km s _1 ) and even broader wings (up to 70% of the total Ha flux, 
with a ~ 290 km s^ 1 ) in regions spatially offset from the clumps by ~ 2 kpc. The broad emission 
likely originates from large-scale outflows with mass outflow rates from individual clumps that are 
l-8x the SFR of the clumps. Based on emission line ratio diagnostics ([Nil] /Ha and [SII] /Ha) and 
photoionization and shock models, we find that the emission from the clumps is due to a combination 
of photoionization from the star-forming regions and shocks generated in the outflowing component, 
with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6xl0 7 -10 8 
cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm~ 3 
in and around the clumps, and ionized gas column densities of 1200-8000 M©/pc 2 ), and SFR (10-40 
M© yr _1 ), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf 
irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the 
nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of 
their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps 
are quite different in terms of their internal properties, energetics and relative ages, and thus we are 
given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant 
clumps at high-z. 

Subject headings: galaxies: high redshift - galaxies: evolution - galaxies: emission lines - galaxies: 
star formation - ISM: jets and outflows 
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1. INTRODUCTION 

Galactic-scale outflows are ubiquitous in lo- 
cal starburst galaxies and in mos t high-redshift 
star-fo r ming galaxies (S F Gs) (iHeckman et al 



199a iLehnert fc Heckmanl 119961 



199a iFranx et all 119971: iPettini et al I 



Steidel et al 
200a 120011: 



Shaplev et all 120031: lAdelberger et~aiT I2003L j200a 



Veilleux et all 12 005: Er b et al. 12006: Wcincr et al. 2009; 
Sato et all 120091 : iRubin et all l2010t lLaw et all 120111 ) . 
These superwinds have a profound impact on the gas 
content and metallicity of the inter-galactic medium 
(IGM), injecting enriched, hot and cold gas into the 
halos and often o ut of the gravitational potentials of 
the host galaxies (iMartinl 119981: iDeyine fc Ballvl 119991: 
Lehnert et all 119991: iStrickland et al l 12004 iRupke et al l 
2005 t iMartinl 120051 : ISteidel et all l20Ta iSturm et all 
2011). Star-formation feedback in the form of these 



winds is key to understanding galaxy evolution, as 
they can mo dulate the star formation activity in 
gas-r ic h disks flKatz et all 119961: iSomerville fc Primac 
1991 iCole et all 120001: iSpringel fc Hernauistl 1200 
' Finlator fc Pavel 120081 : iKeres et all 120091 : iGenel et all 
20101: iDave et^l201ll) . 



Galactic winds in local starburst galaxies are most eas- 
ily observed in edge-on systems in optical, near-IR, soft 
x-ray and radio emission, but outflows have also been 
detected in face-on galaxies from blue-shifted absorption 
by dense cloud s in th e ISM in front of the galactic disks 
(|Veilleux et al.l 120051 and references therein). While 
these winds are clearly important for ongoing galaxy 
evolution and star-formation, little is known about the 
role of such large-scale outflows in high-redshift galax- 
ies, particularly between z ~ 1-3, when massive local 
galaxies were undergoing their most powerful epoch of 



star formation and nuclear activity dMadau et al 



Lilly et al.l [1996] ISteidel et all 119991: lJuneau et al 
Le Floc'h et al.M2005l: iLe Borgne et alll2009D 



1996 



2005 



In the last 15 years, outflows at high redshift have 
been increasingly observed, most often in Lyman Break 
Galaxies (LBGs) and QSO absorption line systems 
(ISteidel et al.lll996l;|Premx^et allll997t IPettini et all | 2000l 
20011: iShapley et all 120031: lAdelberger et all 120031 120051: 
Veilleux et al l 120051: IWeiner et all l2009t faubin et al l 
2010D . Recently, Shapiro et al. (2009) observed broad 



Ha emission lines from the stacked spectra of 47 z ~ 
2 SFGs, possibly indicative of outflows of several hun- 
dred km s _1 and mass outflo w rates comparable to the 
galaxies' SFRs. In addition, ISteidel et all (|2010D found 
interstellar (IS) absorption and Lya emission from gas 
spatially offset from ~100 z ~ 2-3 LBGs using back- 
ground galaxies, indicating that most of their sample 
of high-z galaxies produced metal-enriched galactic out- 
flo ws. A sub s et of this sample was further analyzed 
by lLaw et al.l (|2011[ ) with HST/WFC3 imaging who 
found that outflow characteristics were correlated with 
rest-frame optical morphology. These previous studies 
did not have sufficient spatial resolution to resolve the 
source of the outflows, however. More recently, with 
high spatia l resolution adaptiv e optics integral-field ob- 
servations, iGenzel et al.l (|2011|) have observed powerful 
outflows emanating directly from massive, M ~ 10 9-10 
Mq, star-forming clumps in five rotationally-supported 
SFGs at z ~ 2. These clumps are at the high-mass end 



of the types of clumps observe d with HST imaging by 
" f2011allH) in z ~ 2 SFGs and by 



Forster Schrei ber et al 

Elmegreen et all (12004 I2005D : l Elmegreen fc Elmegreen 
( 2006): lElmegreen et al l (|2009D ; Elmegreen fc Elmegreen 
(2010): iGuo et al.l ([201 ID in tadpole, chain, and clump- 
cluster galaxies in the Hubble UDF at z ~ 1-5. In this 
paper, we present deep, follow-up, multi-band observa - 
tions of one of the galaxies from IGenzel et all (|2011|) . 
ZC406690. 

ZC406690 is part of the "zC-SINF" program aimed 
at studying the spatially-resolved kinematic and star- 
formation properties of z ~ 2 galaxies with the SINFONI 
integra l field spectrograph at the ESO Very Large Tele- 
scope (|Mancini et alll201lD . The galaxy was originally 
draw n from the K-band limi ted sample in the COSMOS 
field ([McCracken et alll2010D with the additional "sBzK" 
color criterion, se lecting candidate 1.4 < z < 2.5 star- 
forming galaxies (jDaddi et al.H2004D . Its redshift of 2.19 
was confirmed from the zCOSMOS-D eep optical spec- 
troscopic campaign (|Lillv et alll2007D . The SINFONI 
seeing-limitcd and adaptive optics-assisted observations 
revealed that it is a rotating, clumpy ring with M* ~4x 
10 ln M Q , and M gQS ~ 2 x 10 11 M B , where M oas is uncer- 
tain b y a factor of 2-3 (|Mancini "etlllll^rGenzel et all 
120111) . With a UV-derived star-formation rate (SFR^y) 
of 337 M yr-\ ZC406690, falls on the SFR-M*main se- 
quence on the high-mass and hi gh-SFR end of the range 
of SINS and zC-SINF galaxies dForster Schreiber et all 
12001 iMancini et al.l l201ll: iWuvts et all 1201 lal) and the 
sample from iForster Schreiber et al.l (|2011al lb1) . This 
SFR is in good agreement with the total amount of star 
formation estimated from Ha corrected for extinction 
with the best-fit Av from SED modeling and adopting 
an extra nebular c orrection of 2.5 (354 M Q yr _1 ) (see 
iCalzetti et aLll2000h. 

IGenzel et all (|201lD found that ZC406690 is ejecting 
mass several times faster than it is forming stars (M out ~ 
3-8 x SFR), as traced through broad, blue-shifted 
Ha emission line wings centered on most clumps, in- 
dicative of an outflowing warm ionized component. This 
finding is consistent with previous galaxy-i ntegrated ob- 
serva t ions of outflows from SFGs at high-z (ISteidel et all 
120101: IPettini et alll2000l: IWeiner et"al]|2009D . With an 
observed wind outflow velocity of 440 km s _1 , the bright- 
est and most massive clump in ZC406690 falls very 
cl osely on the v OM t-M »and v ou t-SFR relations observed 
by IWeiner et al.l ((2009) for high-z galactic outflows. This 
suggests that whatever is driving the relation between 
star-formation and feedback on galactic scales still holds 
on the kpc-scales of giant clumps. 

This compelling evidence for massive, galactic scale 
outflows originating in ~1 kpc clumps has motivated 
us to follow-up on these original observations. We ex- 
ploit here the deep SINFONI+AO data to create indi- 
vidual spectra of not just the clumps, but of the sur- 
rounding regions, whose signal is dominated by the out- 
flowing component projected onto the disk. With ex- 
cellent S/N we are able to create emission line maps 
of [NII]A6584 and [SII]A6716 as well as the broad and 
narrow components of the Ha emission line, enabling 
us to probe the spatially resolved structure and density 
of the wind. With new seeing-limited J and H band 
data of the galaxy, we obtain emission line fluxes of 
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[OIII]A5007,4959, H/3, and [OII]A3726,3729. This new 
data enables us to put the clumps in the BPT d iagram 
(|Baldwin et al.lll981l;IVeilleux fc Osterbrockll987ft . along 
with local SFGs and starbursts as well as other high-z 
SFGs. We can then compare the ISM properties (ion- 
ization state, SFR, density, metallicity, star-formation 
history) of spatially-resolved high-z star-forming clumps 
with those of different types of local star-forming regions 
to better understand the conditions required for these 
massive outflows and how the clumps fit into the overall 
picture of galaxy evolution. 

In section 2, we outline the observations used as well as 
our reduction and analysis techniques. In section 3, we 
present the new spectra of individual regions, emission- 
line maps, and shock-diagnostic diagrams based on more 
detailed examination of the K-band data and the new 
H- and J-band data. In section 4, we compare the re- 
sults to those of local and high-z star forming galaxies 
and discuss the implications. In section 5, we summa- 
rize our findings. Throughout this paper, we adopt a 
ACDM cosmology with Q m = 0.27, fl h = 0.046, an d 
H = 70km s- 1 Mpc" 1 (WM AP7, IKomatsu eTld1 l20Ti"l). 
as well as a IChabrierl (|2003l ) initial stellar mass function 
(IMF). 



2. OBSERVATIONS, DATA REDUCTION AND 
ANALYSIS TECHNIQUES 

2.1. Observations and Data Reduction 

We observed ZC406690 wit h SINFONI/VLT inte- 
gral field unit sp ectroscopy (jEisenhauer et al.l 120031 : 
iBonnet et ai7 2004) using natural guide star adaptive op- 
tics (AO) in K band for 10 hours to achieve a spatial 
resolution of 0.22" (1.85 kpc) FWHM (as presented in 
iGenzel et al.l (|2011| )). Each observation block (OB) con- 
sisted of an "O-S-O-O-S-O" sequence, with 4 on-source 
exposures and 2 others taken on empty 'sky' positions 
about 20 arcsec away from the target. For each on-source 
frame, the object's position was varied around the center 
by ± 0.4" to avoid redundancy. Off-source "sky" frames 
were necessary because the source fills a significant frac- 
tion of the 3.2" x 3.2" FOV at that pixel scale. Individual 
exposure times were of 600s, giving 2400s on-source per 
OB, and thus a total on-source integration time of lOh 
(for 15 OBs). The resulting cube has 0.05" spatial pixels 
and 0.000245 /im spectral pixels. 

We also observed the galaxy in seeing-limited mode 
in J- (100 minutes) and H-band (120 minutes), obtain- 
ing -0.6" FWHM with a pixel size of 0.125" x 0.000195 
[im. With this pixel scale, the larger FOV of 8" x 8" 
allowed an on-source dithering strategy, such that each 
OB consisted of six on-source exposures, with the source 
positioned alternatively in one or the other half of the 
FOV, with additional small dithering of —1/10 of the 
FOV to avoid redundant positions. Individual exposure 
times were all of 600s, giving 3600s on-source per OB, 
and a total on-source integration time of 2h per band 
(for 2 OBs each). Note that we only use 100 minutes of 
integration time in the reduced J-band cube, as the first 
two exposures had faulty pointings. 

For all AO and no- AO OBs, images of the acquisition 
star taken before the science exposures were used to mon- 
itor the positioning and the PSF for the OBs. Standard 
stars of O or early-B types were observed each night, at 



similar airmass as the science data, and used to correct 
for the telluric transmission and to provide the absolute 
flux calibration. 

For th e data reduction, we u sed the software package 
SPRED (jSchreiber et al.l 12004] ) and custom routines for 
optimizing the background/OH airglow subtraction. For 
combining the individual exposures, we determined the 
relative offsets between frames using the known offsets 
for frames within individual OBs and using the position 
of the aquisition star, for frames in different OBs. The 
effective PSF FWHM was determined from the combined 
PSF data sets associated with all OBs in a given band 
(fitting a 2D gaussian profile to the image made by aver- 
aging over all spectral c hannels). For more on the dat a 
reduction procedure, see lForster Schreiber et all (|2009ft : 
iMancini etall (|2011h . 

We calculat e the SFR from the Ha luminosity using 
thelKennicuttl (119981) conversion adjusted to our adopted 
IChabrierl (120031) " IMF, and corrected for extinction. The 
extinction correction is described in Section 3.1. 

2.2. Emission Line Maps 

With the reduced K-band data cube, we create Ha, 
[NII]A6584, and [SII]A6716,6731 maps by simultaneously 
fitting Gaussian components to each of the emission lines, 
holding the nebular [Nil] A6584/6548 ratio fixed at the 
atomic value (3.28) and using the kinematics (line center, 
dispersion) derived from Ha and errors determined from 
the noise cube. Prior to fitting the emission lines, we 
smooth the data by 3 spatial pixels in each direction and 
3 spectral pixels (5x3 for the [SII] map) to boost the S/N 
for the fainter nebular lines. Note that we only include a 
map for [SII]A6716, since [SII]A6731 is strongly affected 
by an OH sky line in some spatial elements. For all line 
maps, we show only pixels with S/N greater than 3er. 

We use the same reduction procedure for the H- and J- 
band cubes, but use a different technique to produce the 
emission- line maps. Since the S/N is not high enough 
to fit each line individually for each spatial pixel, the 
maps were obtained by spectrally integrating over the 
full width of each line. We further increase the S/N of 
the maps by smoothing over 3x3 spatial pixels. After the 
smoothing, all three maps were clipped by a S/N ratio 
of 3 and were overlayed with [OIII] intensity contours. 

Broad and narrow Ha emission line maps are created 
by fitting a two-component Gaussian to the emission 
lines. The resulting broad Ha map is used to visual- 
ize the trend in broad emission surrounding one of the 
clumps, but the map itself has very low S/N, so we gen- 
erate an alternate map for display purposes. This other 
broad Ha map is created by spectrally integrating over 
the wings of the Ha line for the velocity intervals [-430, 
-140 km s^ 1 ] and [280, 570 km s^ 1 ] around the centroid 
of the narrow component. The two maps of the broad 
Ha emission line are qualitatively similar. 

2.3. Extraction of Spectra from Individual Regions 

We create spectra of individual regions from our K- 
band data (Figures 1 and 2) by averaging the spectra over 
both narrow-line (clump) and broad-line (wind) regions. 
We select the wind component based on the region of el- 
evated [Nil] /Ha around the clumps (see section 3.3). We 
fit the broad emission lines using two-component Gaus- 
sian fits with the velocity and width of all lines held equal 
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to that of the Ha line. In some cases, the lines were fit 
equally well with a single Gaussian, which we assumed to 
be the case if the difference in \ 2 values over the Ha and 
[Nil] region of the spectrum between the one- and two- 
component fits was less than 10%. For those spectra 
(clumps C and D), we use the single Gaussian fit. We 
find that the emission line profiles in all regions are well 
described by either one or two Gaussians. We report the 
relative contribution to the flux by the narrow and broad 
components in the next section and Table 1. We derive 
the error for the fit parameters (velocity dispersion, ve- 
locity, etc) by generating 1000 monte carlo realizations 
of the fit and taking the standard deviation for each pa- 
rameter. 

For the H- and J-band spatially integrated spectra, we 
fit each of the [OIII]A4959,5007, H/3, and [0II]A3726,3729 
lines with a single Gaussian component because these 
lines have too low S/N to fit the underlying broad com- 
ponent. For H-band, we fit all of the lines using the 
kinematics of the [OIII]A5007 line, and for J-band we 
constrain the two [Oil] lines to have the same kine- 
matics as Ha and fix their wavelength separation to the 
known value, but allow their relative intensities to vary. 
For emission line ratios that contain lines from multi- 
ple bands, we employ an extinction correction as out- 
lined in Section 3.1 with an E(B-V) as determined from 
Ha/H/3 for clump A. In addition, for all ratios (except 
[Nil] /Ha and [Sll]/Ha), we use emission line fluxes de- 
rived from single-component Gaussian fits. When com- 
paring K-band (higher spatial resolution) and H-band 
emission lines, we resample the K-band cube to match 
the H-band cube, and convolve the new K-band cube 
with the H-band PSF. The H- and J-band cubes have 
the same spatial and spectral resolution. We therefore 
use the exact same spatial pixels when extracting spec- 
tra of individual regions from H, K and J bands. 

To create the [SII] line profiles (Figure 2) and ratios, 
we generate a velocity-shifted cube to increase the S/N. 
This is done by shifting the spectrum of each pixel based 
on the measured Ha velocity such that all pixels have 
their Ha peak at the same wavelength. This method 
has bee n shown to increase the S/N of region-integrated 
spectra (|Shapiro et al.l l2009) . To negate the effects of an 
OH sky line that lies near the [SII]A6731 emission line 
for parts of the galaxy, we create clump A and clump B 
spectra that have maximal S/N ratios in the wavelength 
range of this line. To this end, we define the spatial 
regions of the clumps by only adding spatial pixels to 
the spectra such that each additional pixel increases the 
S/N over a given threshold. Essentially, we only add 
pixels with a [SII]A6731 S/N ratio > 4. In this way, we 
only use measurements of the [SII]A6731 line in pixels 
that are minimally affected by the OH sky line. In or- 
der to maximize the S/N, these regions must be quite 
extended, so to get [SII]A6731 measurements for smaller 
areas within these larger regions (i.e. for the broad re- 
gions surrounding the clumps), we assume the measured 
[SII]A6716/[SII]A6731 ratio within the larger regions is 
constant. We can then estimate [SII]A6731 fluxes based 
on the measured [SII]A6716 flux in the smaller regions 
and the measured [SII] ratio from the corresponding 
larger region. 

The H/3 line is obstructed by an OH atmospheric line 
for about half of the galaxy. Therefore, to obtain an 



H/3 flux estimate for clump B (which falls in this right 
half), we assume clump B follows the extinction of the 
dense, clump gas and has the same observed Ha/H/3 ratio 
as clump A. We obtain an H/3 flux using this ratio and 
the Ha measurement for clump B, although this estimate 
is much more uncertain than the clump A H/3 flux. 

3. RESULTS 

3.1. The Clumps are More Highly Extincted than the 
Surrounding Gas 

The extinction correction used to calculate the SFR is 
based on the best-fit visual extinction derived from stel- 
lar population modeling of the UV to mid-IR spectral 
energy distri bution (SEP) presented by iMancini et al.1 
(IMl L for alCalzetti et al l (pOOOT ) reddening law. Fol- 
lowing rCalzetti et al.l ()200*oT ). we assume that the Ha line 
emission is more attenuated than the bulk of stel- 
lar continuum light and thus correct the observed 
Ha luminosity for Ayi^htd^ = A-viSF.m x 2.5 (see 
also iForster Schreiber et al] [20091 ; IMancini et"aT1 120111 : 
iWuvts et al.H2011al ). 

We also calculate the extinction for clump A and the 
left-half of the galaxy (which is unaffected by an OH sky 
line ne ar the H/3 line) based on the Ha /H/3 ratio and us- 
ing the lCalzetti et all (12000T) reddening l aw, assuming an 
intrinsic ratio of 2.86 (jOsterbrocldll989l ) appropriate for 
an electron density of 100 cm~ 3 and a temperature of 10 4 
K. We fin d that the SED-base d Ay (described above) as - 
sumed in iGenzel et all (|2011[ ) and IMancini et~all (|2011l ) 
of 1.3 ± 0.5 (E(B-V) = 0.33 ± 0.1) is consistent with 
the diffuse emission from the galaxy (as calculated with 
Ha/H/3 from the left-half, A v = 1.5 ± 0.8, E(B-V) = 0.37 
± 0.18) and that the nebular emission in clump A is more 
highly extincted (Ay = 2.0 ± 0.8, E(B-V) = 0.49 ± 0.17), 
as one would expect for a star-forming region. All of the 
forthcoming calculations use the clu mp Ha flux a s calcu - 
lated and corrected for extinction in lGenzel et al.l (|2011[) . 
If we were to apply an extinction correction based on 
the Ha/H/3 measurement for clump A, the main clump 

properties (Lna, SFR, M out , S gQ s) would be increased 
by about 50%, and therefore the calculations made in 
Section 4 based on these quantities would be scaled ac- 
cordingly. This relatively small factor will not affect any 
of the conclusions made in this paper. 

3.2. Spectra of the Clumps and their Surroundings are 
Indicative of Localized Superwinds 

We explore the properties of the star-forming clumps 
by creating individual spectra of both the clumps and 
their surrounding regions. In the regions surrounding 
the clumps, we find broader Ha emission lines and en- 
hanced [Nil] /Ha and [SII] /Ha ratios as compared to the 
clump regions. Enhancement of these emission line ratios 
could be due to several physical processes including the 
presence of shocks in a large-scale outflow or an AGN. 

In this case, however, it is very unlikely that these 
large line ratios are due to an AGN. First, the region in 
which we find these large [Nil] /Ha and [SII] /Ha ratios 
is extended, unlike an AGN point source. Second, these 
regions are located around clumps in the disk (or ring) of 
the galaxy, far from the center, where one would expect 
to find an AGN. The velocity field of the ring is also very 
regular and characteristic of face-on rotation, unlike what 
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Fig. 1. — Panels a, b and c show K, H and J band spectra of clump 
A with the relevant emission lines labeled, respectively. The black 
line represents the data, while the red line shows our best fit. The 
wavelength axis is in rest-frame units. 

would be expected if clump A and clump B are the nuclei 
of merging galaxies. Third, the value of the elevated 
[Nil] /Ha ratio in this region (0.3-0.5) is much less than 
those expected for AGN (at least 0.6-0.8). Fourth, we 
find that the [Nil] emission is well fit by both a broad 
and a narrow component (see later in this section), which 
is inconsistent with classical broad-line regions around 
AGN. Finally, there is no evidence of an AGN from the 
ZCO SMOS optical spectra of this galaxy ()Mancini et al.l 
USB- 

We further discuss the possibility of shocks in Section 
3.3, but even without the emission line ratios, the broad 
emission line wings (up to 800 km s _1 ) in the spectra of 
the regions surrounding the clumps seem to indicate the 
presence of outflowing material. We therefore call these 
regions the clump A and clump B 'wind' regions. 

In Figure 1, we show K-, H-, and J-band spectra from 
Clump A (the brightest and most massive clump) along 
with our best fit model. In Figure 2, we show emission 
line profiles of Ha and [NII]A6584 for the four clumps 
in ZC406690 and the wind regions surrounding clumps 
A and B. Clumps A and B and the surrounding regions 



are well fit by two Gaussian profiles, with narrow compo- 
nents with FWHM = 200±8 km s _1 and FWHM = 210±9 
km s _1 and underlying broad components with FWHM 
= 470±28 km s^and FWHM = 630±82 km s"\ which 
we ascribe to outflowing material. On the other hand, 
clumps C and D can be well fit by one component with 
FWHM = 92±2 km s^and 89±2 km s _1 , suggesting lit- 
tle or no outflow from these clumps. The clump B wind 
spectrum is the most extreme example of a localized su- 
perwind with a very broad blue line wing, with the edge 
of the wing extending out to roughly 800 km s _1 from 
line center. In addition, the centroid of the broad com- 
ponent is blueshifted by 150±16 km s^from the narrow. 
This is consistent with the signature of a powerful out- 
flow, as one would expect to see the majority of the out- 
flowing material as blueshifted if the clumps are heavily 
obscured. The redshifted emission would fall behind the 
plane of the galaxy, and much of this emission would be 
extincted as it passes through the clumps as suggested 
by our Ha/H/3 measurements. 

The fraction of Ha emission that comes from a broad- 
velocity component varies substantially for different re- 
gions in the galaxy. For the galaxy as a whole, broad 
emission accounts for roughly 33% ± 2% of the total 
Ha flux. For clumps A and B and wind A, the broad 
component accounts for about 50% ± 6% of the total 
Ha flux, while for wind B, broad emission comprises 70% 
± 3% of the Ha flux. In all of the clump and wind 
regions, the narrow component likely derives from star- 
formation in the disk, while the broad component may 
come from a combination of turbulence and outflowing 
gas. In the case of wind B, the broad velocity-component 
peak wavelength is blueshifted slightly from the narrow- 
velocity component peak, which indicates that a large 
fraction of the Ha flux from this region comes from an 
outflowing wind. The emission line profile for [OIII] is 
very similar to that for Ha, such that the clump A [OIII] 
line is slightly less broad than Ha, but for clump B, both 
lines are equally broad, with the asymmetric blue line 
wing also present for [OIII]. However, we note that the 
aforementioned comparisons are done by eye since the 
S/N for the H-band spectrum for each clump is too low 
to perform a rigorous two-component fit. Since [OIII] 
and [Oil] likely have similar kinematics, our choice of 
using Ha kinematics to fit the [Oil] lines is thus sup- 
ported. 

We derive the local electron density of the clumps us- 
ing the ratios of the [SII] and [Oil] doublets. The lower 
right and left panels of Figure 2 compare the spectra and 
one-component Gaussian fits of the [SII] emission lines 
for clumps A and B, with regions selected by the method 
described in 2.3. For clump A, [SII] A6716 / [SII] A6731 = 
0.7 +/- 0.1, while for clump B, = 1.2 +/- 0.3, corre- 
sponding to electron densities of 1800± 1000 and 290±300 
cm~ 3 , respectively (|Osterbrocklll989D . These densities 
likely come from a combination of both the star-forming 
regions in the disk and high-density shocked filaments 
or clouds in the outflow, from which the emission lines 
originate. Since the clump+wind B region has a larger 
broad Ha flux fraction than the clump+wind A region, 
we would expect that more of the [SII] emission for 
clump B derives from the outflow than for clump A. 
Therefore, the observed clump A electron density of 1800 
cm~ 3 is probably a better estimate for the dense clump 
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Fig. 2. — Emission line profiles and best fit to Ha and [NIIJA6584 features for six regions of ZC406690: clump A, clump A 'wind', clump 
B, clump B 'wind' and clumps C and D. Included with these spectra are one- or two-component Gaussian fits. For the two-component 
fits, the narrow component is shown in red and the broad in green. For the spectra with one-component fits, these are shown in red. The 
errors (la rms) are shown below in blue. Clumps A and B and their winds are well fit by two Gaussians, while clumps C and D can be 
well fit by one component. The lower right and left panels of Figure 2 compare the spectra and one-component Gaussian fits of the [SII] 
emission lines for clumps A and B, with regions selected by the method described in 2.3. For clump A, [SII] A6716/[SII] A6731 = 0.7 +/- 
0.1, w hile for clump B, = 1.2 +/- 0.3, corresponding to electron densities of 1800 ± 1000 and 290 ± 300 cm -3 , respectively l|Osterbrockl 
1989). The grey hatched lines in all figures show the wavelength range of strong sky OH emission features. The central panel is an Ho map 
of ZC406690, with each region outlined. 



environment, while the clump B density of 290 cm' 3 is 
likely an average of clump gas and less dense outflowing 
cloud gas. The clump B estimat e is similar to the value 
assumed bv lKewlev et al.1 ()2001al ) in their starburst mod- 
els (350 cm -3 ), which is based on the average density of 
individual HII regions from their sample observed with 
a 1 kpc slit. 

Densities calculated from the Oil doublet are much 
more uncertain (as the lines are mostly blended together 
in our spectra), and are only consistent (given the un- 
certainty) with the [SII] derived densities for clump B 
(see Table 2). This discrepancy could could be due to the 
two tracers probing different regions in the wind. The la 
ranges quoted in Table 2 for both the [SII]- and [OII]- 
derived densities are obtained by adding and subtracting 



the uncertainties derived from the emission line fits to 
the ratio valu e and finding the corresponding electron 
densities from lOsterbr ock (1989|). 

We interpret the wind B region as the projection of an 
outflowing component onto a low-inclination disk (~ 20° 
(ts ), as estimated in iGenzel et~aT1 (|2011| )). bmce re- 
gions adjacent to the clumps would have much lower 
Ha luminosity from star formation than the clumps, the 
broad emission from the winds can account for much 
more of the total signal. Furthe rmore, if we assume the 
outflow has a biconic geometry (jVeilleux et alJl2~005L and 
references therein), then the distant parts of the out- 
flow will have a larger spatial extent than the clump, 
further separating this "wind" region from the clump in 
projection. This notion of viewing the wind in projec- 
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tion is consistent with the presence of less of an offset 
broad component for clump A, as it is likely younger 
than clump B (see section 4.2) and therefore has had 
less time to produce an extended wind. The relative 
outflow velocities from the two clumps further supports 
this view (with \ mi *B >> Vnut.A), as most ou tflow mod- 
els (jHeckman et alll990l:IVeilleux et alj|2005| ) predict in- 
creasing wind velocity with distance from the galactic 
plane. 

We can use this model of a projected wind to learn 
about the geometry and structure of the outflow. The 
offset between clump B and wind B is ~ 1.8 kpc, so 
with an inclination of 20°, the average distance of the 
line-emitting outflow from the disk is ~ 5.3 kpc. Given 
that the broad component of wind B is blueshifted by 
~ 175±20 km s _1 from the narrow component of clump 
B, then the average timescale of the line-emitting part of 
the wind is 30 Myr. It is likely that the gas in clump B 
has been forming stars for much longer than this given 
its average metallicity (see Table 2), but we are only 
observing the portion of the outflow coming from dense, 
hot filaments in the wind, and the older parts of the 
outflow may be too diffuse or faint to observe. 

3.3. Broad Emission is Spatially Correlated with 
Shock-like Emission Line Ratios 

Figure 3 shows spatially-resolved emission line maps 
for Ha, [NII]A6584 and [SII]A6716, which illustrate the 
effect of the outflowing gas on the region surrounding 
the clumps. The peak in [NII]A6584 emission is located 
near the peak of Ha emission over the clumps but the 
[Nil] emission is more extended around clump B than 
the Ha emission (by ~1 kpc). This effect is seen even 
better in the ratio map of [Nil] /Ha located in panel e of 
Figure 3. There is a semi-circle like structure centered on 
clump B that peaks in [Nil] /Ha over the broad-line re- 
gion, with an enhancement of about 80% from the clump 
value (see Table 1). The same effect is seen around clump 
A with a smaller enhancement. Errors for the different 
regions are over an order of magnitude lower, such that 
the uncertainty ~ 2-5% on and near the clumps. These 
enhanced values of the [Nil] /Ha ratio can be indicative 
of both hi gher metallicity and eviden ce of shocks and 
AGN (e.g. iDopita fc Sutherland! I1995Q . We previously 
discussed why an AGN explanation is unlikely in Section 
3.2. 

We can somewhat disentangle the effect of metallicity 
and shocks on the [Nil] /Ha ratio by examining other 
emission line ratios and comparing the data to shock 
and photoionization models from the literature. Both 
the [SII]A6716+6731/Ha and [OlII]/H/3 ratios are much 
more sensitive to the ionization parameter (and thus 
shocks) than to metallicity, in the metallicity regime of 
the clumps. 

The [SII]A6716 map (panel c) shows that much of the 
emission is located on or near the clumps, although it 
does not follow the Ha emission very closely. Panel f 
shows a map of [SII]A6716/Ha, which serves as a lower 
limit for [SII]A6716+6731/Ha which is omitted due to 
the aforementioned OH sky lines. Even with this un- 
derestimate, we see slightly higher ratios surrounding 
the clumps than over their centers (0.059 vs. 0.087 for 
clump A and 0.11 vs. 0.14 for clump B). This is consis- 
tent with the notion that the larger [Nil] /Ha ratios are 



caused by shocks (and possibly also photoionization from 
star formation) as opposed to sharp metallicity gradients. 
The typical uncertainty for the [SII]A6716/Ha ratio near 
the clumps is 5-10%, although the S/N drops off sharply 
away from the clumps. 

It is possible that there is some variation in the 
[SII]A6716/6731 ratio between the clump and wind re- 
gions that we are unable to measure because of the 
OH sky lines, and we have considered the scenario in 
which this variation could compensate for the gradient of 
[SII]A6716/Ha from the clump to wind regions, thereby 
erasing the overall [Sll]/Ha ratio gradient. However, if 
we calculate the variation in the [SII] ratio required to 
erase the [SII] /Ha gradient and assume the clump re- 
gion has the observed value of the ratio, we get implau- 
sible values for the [SII] ratio of the winds, 1.6 for the 
clump A wind and 2.5 for the clump B wind, which are 
inconsistent with the co l lision al de-excitation model at 
10 4 K from lOsterbrockl (| 19891 ) . Thus the gradient in 
[SII]A6716/Ha between the clump and wind regions im- 
plies a similar gradient in [SII]A6716+6731/Ha, which 
cannot be washed out by density variations. 

This region of enhanced emission line ratios surround- 
ing clump B is spatially correlated with the bulk of the 
broad Ha emission, which we ascribe to an outflowing 
component. From the broad Ha map (Figure 3 panel d), 
it is evident that most of the broad emission comes from 
directly over clump A and to the SE of clump B, coin- 
ciding with the enhanced [Nil] /Ha and [SII] /Ha ratios. 
This broad region near clump B, from which we extracted 
the 'wind B' spectrum, not surprisingly also shows up 
as a peak in the velocity dispersion in panel g, and is 
blueshifted relative to clump B, as seen in the velocity 
map (panel h). The latter statement is confirmed by ex- 
amining the line centers of the narrow components for 
clump B and wind B. The kinematic maps are derived 
by fitting a single G aussian comp o nent to the Ha line 
(for more on this see iGenzel et all (|2011[ )). Thus, if ac- 
tive star formation in the clumps produces a large-scale, 
energetic outflow, this could explain both the broad emis- 
sion lines and the enhanced emission line ratios in this 
region. 

In order to visualize the emission line gradients around 
clump B and understand their errors, we show the aver- 
age flux across a swath of the emission line maps with 
a slit width of 0.2" in Figure 4. This figure shows that 
there is a relative increase in emission for both [Nil], 
[SII] and Ha broad relative to Ha narrow around clump 
B with errors from the noise cube that are much smaller 
than these gradients, and thus the emission line gradients 
seen in Figure 3 near clump B are real. While the errors 
increase dramatically for [SII] /Ha near the edge of Fig- 
ure 4, the gradient is seen clearly near clump B where the 
S/N is relatively high. The Ha broad emission is calcu- 
lated as described in Section 2.2 using a two-component 
fit to the Ha line. 

In Figure 5, we show the emission line maps for [OIII], 
H/3, and [Oil]. As is clear from the first and last panels, 
the [OIII] and [Oil] emission is coincident with the lo- 
cation of the clumps. Even with the final smoothing (see 
Section 2.2), the S/N in individual pixels for H/3 is very 
low, and thus only a small region near clump A remains 
after the S/N = 3 clipping. Fortunately, the H/3 line in 
the integrated spectrum for clump A has a S/N ratio ~ 
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Fig. 3. — Spatially-resolved emission line maps for (a) log[H« narrow], (b) log[[NlI]A6584], (c) log[[SlI]A6716], (d) log[Ha broad], (e) 
log[[NH]/Ha], (f) log[[SlI]A6716/Ha], (g) the Ha dispersion and (h) the Ha velocity field. The top three panels show the velocity-integrated 
flux and have a S/N per pixel (for the un-resampled maps with 0.05" pixels) near the clumps of 100-200, 10-40, and 10-40, respectively. 
[SIIJA6731 is not shown since it is affected by an OH sky line. All of the maps are rebinned by two pixels and smoothed by 3 pixels for the 
purpose of presentation and are overlaid with Ha intensity contours to highlight the locations of clumps A, B, C and D (labeled in panel 
a). The rebinning does not change the overall appearance of the maps as the S/N is very high, but merely accentuates features that are 
already present before rebinning. The spatial resolution of these maps is 0.2". From panels (d), (e), and (f), we find that broad emission 
surrounding clump B is spatially correlated with shock-like [Nil] /Ha and [Sll]/Ha ratios, indicating an energetic wind coming from the 
clump. 
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Fig. 4. — Averaged flux of several emission lines and ratios along 
a four pixel- wide (0.2") swath of the galaxy moving from clump D 
to B to C. (a) Average flux of Ha narrow, Ha broad, [NIIJA6584 
and [SII]A6716, along with error bars from the noise cube, (b) 
[Nil] /Ha and [SII] A6716/Ha and error bars from the noise cube 
with Ha narrow scaled in the background for reference. There is a 
relative increase in emission for both [Nil] , [SII] and Ha broad rel- 
ative to Ha narrow around clump B and the errors are much smaller 
than these gradients. For [Nll]/Ha, [Sll]/Ha and Ha broad, the 
gradients from the peak of emission (or the peak of the ratio) to 
the minimum (near the clump B center) are over 10a. The inset in 
the upper left corner of (a) shows the width, location and direction 
of the swath. 



5, and we are able to use this to derive [OIII] /H/3 ratios 
for the clumps. 

Fig ure 6 shows the BPT diagram (Baldwin et al. 
I1981D with SPSS data from the latest DR7 release 
(jAdelman-McCarthv fc et all [20091) along with the data 



for clumps A and B, to compare the latter to local star- 
forming galaxies and other high-z galaxies. Since the 
clump A and B data points use both high-resolution K- 
and low-resolution H-band data, we extract larger re- 
gions, which contain both the 'clump' and 'wind' regions. 
We are unable to get a direct measurement of the H/3 line 
for clump B, as it is shifted onto an OH sky line in that 
part of the galaxy, but we estimate the clump B H/3 value 
assuming a constant Ha/H/3 ratio for the clump regions 
(Ha/H/3 = 5.1 as measured from clump A). The position 
of the clump B point is therefore much more uncertain, 
and the lower error bar for [OIII] /H/3 extends off the fig- 
ure. The errors come from the uncertainty of the fits to 
the emission lines. Clumps A and B fall in the composite 
region of the BPT diagram that is generally occupied by 
starburst and high-z galaxies, which are generally found 




Fig. 5. — Spatially-resolved emission line maps for (top to bot- 
tom) [OIII], H/3, and [Oil]. The maps have been smoothed by 3x3 
pixels to enhance the S/N and then clipped by S/N = 3. The maps 
were then rebinned by 2 spatial pixels and smoothed again by 3 
pixels for the purpose of presentation. The white contours trace 
the [OIII] flux. The un-resampled cubes have 0.125" pixels and a 
PSF FWHM of 0.6". 

to have higher SFRs, densities and ioniz ation parame- 



ters than normal SFGs and HII regi ons ( Shapl ev et al 



| 2005| lErb et al.ll2006l : iLiu et al.l f2008: Br inchmann et al 
120081) . Their offset could also be due to shocks, AGN, 
or a different N/O ratio. As the clumps occupy the 
same parameter space as these galaxies, it appears that 
their star-forming regions more closely resemble those of 
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Fig. 6. — The BPT diagram including data from several high-z 
galaxies along with the clumps from ZC406690. Clumps A and 
B are shown as a red asterisk and blue square, respectively. The 
lower error bar for clump B extends beyond the figur e boundaries . 
Other hig h-z data is sh own as grey squ ares (Shaplcy ct al, I2005f ), 
circles l|Erb et alJI^OrD . and diamonds QLiu et al.1120081 ). A subset 
of SDSS galaxies are shown in g rey with the empi r ical A GN/HH 
boundary as the dashed line from Kauffmann ct al. (2003) and the 
theoretical boundary as the solid line from Kcwlcv ct al. (2001b). 
The clump data points fall in the region offset from the main HII 
branch that is also occupied by other high-z galaxies and many 
local starbursts (as shown by ILiu eF al. 2008; Brinchman n et all 
120071 ). 



starbursting galaxies than of HII regions in normal local 
SFGs. 

In order to differentiate the effects of metallicity, ion- 
ization parameter, and the contribution of shocks, we 
compare our data to phq toionization and shock mod- 
els f rom IRich et al.l ( 2011 ) and photoionization models 
from lLevesque et al.l ( 2010 ). These photoionization mod- 
els are similar in many respects: they are both cre- 
ated with the code s Starburst99 (iLeitherer et al.l 11999ft 
and Mappings III (jSutherland fc Dopital 11993ft . use the 
same evolutionary tracks and both assume continuous 
star-formation histories, but they also differ somewhat 
in terms of model atmosphe re, geometry and IMF. Al- 
though the IRich et ail (|2011| ) shock models are for slow 
shocks (up to 200 km s" 1 ), they sample a large portion 
of the parameter space (in terms of metallicity and ion- 
ization parameter) which fits our observations, and the 
exact velocity of the shocks do not qualitatively effect 
our conclusions since there is a much bigger difference 
between photoionization and shock models than slow vs. 
fast shocks in terms of emission line ratios. While ev- 
idence for shocks indicates the presence of an outflow, 
emission from both shocks and photoionization could be 
present in a wind. 

Figure 7 shows our data from the clump and wind 
regions as well as the photoionization and shock 
models in the [Nil] /Ha vs. [Sll]/Ha! plane (with 
[SllJ/Ha calculated as described in section 2.3), with 
varying metallicity, ionization parameter and shock ve- 
locity. We are able to get data for both the clump 
and wind regions here as we are only using the high- 
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Fig. 7.— The [Nll]/Ha vs. [Sll]/Ha plane with data from the 
clump and wind regions as well as the photoionization and shock 
models (with [Sll]/Ha calculated as described in section 2.3). The 
data are shown as a red asterisk, a purple square, a blue 'x' and 
a cyan star with error bars (for clump A, w ind A, clump B and 
wind B, respectively). The Rich ct al. (20TJ1) and ILevesque et al.l 
(2010) photoionization and shock models are shown a s the black 
and g rey lines. The darkest /dashed lines are the IRich et alj 
(I201l|) photoionizat i on mo dels (RP), the medium/dotted lines are 
the ILev esque ct al. (2010) phot oionization models (LP) and the 
lightest /dash-dot lines are the IRich ct al. (201J]) shock models 
(RS). Open circles denote the models with O/H = 8.69 and open 
diamonds denote models with O/H = 8.39. For the photoioniza- 
tion models the points along the line represent a varying ionization 
parameter, increasing to the left (logU = 7—8). For the shock-only 
models, the points along the line vary in shock velocity, increasing 
upward (v = 100-200 km s — 1 ). The solid black lines indicate the 
ionization parameter range for each metallicity line as calculated 
in section 3.3. These measured ionization parameters are offset 
from those implied by the photoionization-only models, indicating 
a shock contribution to the emission. 

resolution K-band data while we could only get overall 
clump regions with lower resolution for Figure 6 since 
we used lower-resolution H-band data in addition to K- 
band. Thus the clump regions in Figures 6 and 7 cover 
different spatial areas and therefore they have emission 
line ratios with slightly different values. We also note 
that, as for Figure 6, the error bars come from the un- 
certainty in the emission line fits and are therefore larger 
than the error bars shown in Figure 4, which are derived 
from the noise cube. 

The two photoionization models sh own in Figure 7 
(jRich et all 120 lit ILevesque etafl 120101 ) are roughly the 
same for [O/H] = 8.39, however, the [O/H] = 8.69 line is 
offset to much larger v alues of the [Nil] /Ha ratio for the 
ILevesque et al.l (|2010ft model. This is likely due to one 
of the differences between the models mentioned above. 
At first glance, the data points seem to be well fit by the 
photoionization-only models. Both wind regions appear 
to fall near the same metallicity line as their respective 
clumps, but offset to higher emission line ratios, indi- 
cating a lower ionization parameter. However, we can 
independently calculate the ionization parameter for the 
clumps, and compare that to the ionization parameter 
implied by the photoionization models. 

We calculate the ionization parameter for the clumps 
using two different methods: 1) Using the O32 diag- 
no stic with metallici t y-dep endent fitting curves based 
on iKewlev fc Dopital J2002|) and 2) using the observed 
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Ha luminosity and calculated electron density as follows: 



where Q = \jLyc/^Lyc is the rate of ionizing p hotons, 
Liyc = 16xLfj a (as given by IQsterbrockl [1989) . r is 
taken as the clump HWHM, and n e is the electron den- 
sity calculated from the [SII] doublet. The results are 
shown in Table 2 and the range of ionization parame- 
ters (for the corresponding metallicities) , is denoted by 
the thick black lines in Figure 7. These ionization pa- 
rameters are 3-10 times higher tha n those predicted from 
the photoioniza t ion-on ly models of lRich et al.l (|2011| ) and 
iLevesque et al.l (|2010f ). The data points appear to be 
offset towards the shock-only models, and are perhaps 
better fit by a model with contributions to the emis- 
sion from both photoionization and shocks. In fact, if we 
compare the [Nil] /Ha and [SII] /Ha ratios for the clumps 
along with the corresponding ion ization parameters with 
Figure 10 from lRich et al.1 (|2011l ). which shows the pho- 
toionization + shock models, we find that the clumps 
and their winds are best fit by models with some shock 
contribution. The clumps are well fit by models with 5 ± 
5% and lOj^ % shock contribution to the nebular emis- 
sion (for A and B), and the wind regions are well fit by 
models with 15j^lj% and 30^2q% shocks, where the er- 
ror bars come from fitting the la values of the ionization 
parameter and emission line ratios to the models. 

Based on the offsets from the observed ionization pa- 
rameters, it is reasonable to assume that the elevated 
emission line ratios in the wind regions compared to the 
clumps are due to a larger contribution to the flux by 
shocks, and that the offset of the clumps in the BPT 
diagram from normal SFGs is due in part to shocks gen- 
erated by star formation feedback. The highly turbu- 
lent, fast outflows emanating from the clumps could eas- 
ily produce shocks as the hot superwind fluid collides and 
mixes with cooler ISM gas. 

4. DISCUSSION 
4.1. Properties of the clumps and their outflows 

Based on the observed properties of the clump A 
and B winds, we can perform simple calculations to 
better understand the structure of these outflows and 
their relationship to local and high-z galactic-scale su- 
perwinds. For the following calculations, we use the 
S FRs, gas mass e s and mass outflow rates as calculated 
in IGenzel et all (|2011l ) The SFR was calculated from 
the Ha flux using the iKenmcuttl (|1998l) calibration ad- 
justed for a Chabrier IMF with the extinction correction 
described in Section 3.1. 

The gas masses were derived from the SFRs using 
the Sc hmidt-Kennicutt relation as given in IGenzel et al.l 
(f2010h . There are systematic errors that arise from this 
calculation due to (1) scatter in the intrinsic relation, (2) 
the varying slope and zero point of the relation as quoted 
by different authors, (3) the extrapolation of the relation 
from low-z to high-z, and (4) the application of the rela- 
tion on ^kpc scales. We address each of these issues in 
Appendix A, and ultimately derive an overall systematic 
uncertainty for the gas masses of the clumps of +0.3/-0.8 
dex. 



The mass outflow rates were calculated assuming an 
outflow of warm ionized gas into solid angle ft with ra- 
dially constant out flow rate and velocity as derived in 
IGenzel et al.l (|2011| ) and are quoted in Table 3. Details 
of the model are discussed in Appendix B. We note that 
for our outflow model, we only calculate the mass and 
mass outflow rate of the warm ionized component of the 
wind, and thus our numbers are lower limits. 

4.1.1. The Winds from Individual Clumps Are Very 
Massive and Energetic 

From our mass outflow rate derivation, we also obtain 
the mass of the warm ionized component in the wind 
(M w ) for each clump. We estimate the energy of the 
wind as E w = \ M^x v out 2 , where v out is the maxi- 
mum outflow velocity determined from the width of the 
Ha emission line (see Table 3). Using an average wind 
velocity instead would only change E w by a factor of 2 
and will not affect any of the forthcoming conclusions. 

We find that the mass and energy in the outflows of 
individual clumps lies in the range l-2xl0 8 M Q and 3- 
7xl0 56 ergs, several times greater than what is observed 
in the warm component from entire ULIRG and dwarf 
starburst galaxies (M»,~ 10 5 " 7 M Q and E m ~ 1Q53-55 
ergs) based on numbers compiled by iVeilleux et al.l 
(2005). We caution that our measurements must be 
carefully compared to those of local starbursts, as the 
mass estimates are strongly dependent on the extinc- 
tion corrections used and the quoted local starburst en- 
ergy estimates are taken from dwarf starbursts and spiral 
galaxies, which likely produce less energetic winds than 
ULIRGs. 

The main point that should be gleaned from this sec- 
tion is that these clump-scale outflows are massive and 
very energetic even by local starburst standards. 

4.1.2. Are the Outflows Energy or Momentum Driven? 

Galactic-scale outflows are considered to be driven ei- 
ther by the energy rele ased from supernovae and the UV 
lumin osity from stars (|McKee fc Ostrik"erlll977l ; iMatznerl 
2002) or through momentum transport by radiation pres- 
sure on dust a nd the momentum from stellar winds 
and supernovae (|Murrav et al.ll2005l [201fl . Momentum- 
driven winds tend to dominate in dense regions where 
much of the energy is radiated away, but the momen- 
tum remains (|Murrav et al.ll2005D . IMurrav et all (j2005T ) 
make several predictions for what one would expect 
given energy- or momentum-driven outflows from star- 
formation feedback in terms of M out (the mass outflow 
rate), p ou t(the momentum outflow rate) and rj (the mass 
loading). 

For a starburst in whic h 90% of the s u perno vae (SN) 
energy is radiated away, IMurrav et al.1 (|2005l ) predict 
that EsN,out ~ 10 _3 Lh o /, as 1% of the total luminosity 
released throughout a stars lifetime is ejected in the 
supernova, and 10% of the supernova energy is not 
radiated away. Thus, assuming the supernovae energy 
goes into the wind, E w j n d/L;, / ~ 10 -3 . For the clumps, 
we find that 

E w ind/L bo ; ~ 5 M out x v out 2 / L fco/ = 0.005 and 0.05 
with ranges of 7.6xl0" 5 -0.05 and 0.0002-0.3 for clumps 
A and B respectively. 
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Thus for a feedback efficiency of 10%, the clump 
A wind could just barely be explained by an energy- 
driven model. However, the clump B wind is inconsistent 
with a purely energy-driven wind unless none of the SN 
energy is radiated away. If however, the extinction in 
clump B is much greater than for clump A, from the 
disk would be underestimated more than M out , which is 
coming from gas that is farther from the disk. This could 
lower E w i„d/Lf, ; for clump B, bringing it closer to a 
number that could be explained by energy-driven winds. 
In addition, the large error bars for these quantities 
(mostly derived from the large error in M out ) do make 
them roughly consistent with the model, although due 
to the fact that M out is likely an underestimate (see 
Section 4.1), the results are more likely skewed to larger 
values of E wind /L 6o; . 

For a purely momentum-driven wind, IMurrav et al.l 
(2005) predict that p w i n d/Prad ~ 1, and we find that 

Pwind/Prad ~ M out x v out / (L boi /c) = 6 and 30 
with ranges of 0.5-65 and 2-245 for clumps A and B 
respectively. 

This result indicates that clump A can also be ex- 
plained by a momentum-driven wind, but clump B 
cannot (pwind/Prad = 30), although again, the large 
error bars make this conclusion somewhat uncertain. 

There are many factors that coul d be biasing this mea - 
surement. First, as mentioned in iGenzel et al.l (| 2 1 If ) . 
our outflow velocities could be too high by a factor of 
2, depending on the definition chosen for the maximum 
outflow velocity. Second, in an optically thick medium, 
the total effect of radiation pressure could be greater by 
a factor of (I+t/^) due to multiple scattering of infrared 
photons, although it will likely not exceed a factor of 
a few. These two effects will bring down the value of 
Pwind/Prad- On the other hand, our mass outflow rates 
could be underestimated, both as a result of only includ- 
ing the warm ionized component of the wind and also 
due to overestimated wind densities (see next section), 
and this would mean our momentum wind efficiency is 
too low. These uncertainties are represented in the large 
systematic errors for pwind/Prad- 

Our results suggest that the clump A wind could be 
either energy- or momentum-driven but that the clump 
B wind is harder to explain by either driver. However, it 
is possible that current models do not fully encapsulate 
the physics in these extreme high-z star-forming regions. 
In the next section we further compare the relative con- 
tributio ns to feedback by S N and radiation pressure us- 
ing the lOstriker fc Shettvl (|201lD model. We note that 
IMurrav et al.l (|2005f ) predict that' in both energy and mo- 
mentum driven outflow models, the mass loading factor 
(?7 = M out /SFR) ~ 1 and this is roughly consistent with 
our o bservations (77 = 3-8). Similarly, iHopkins et al.l 
(|20TlT) predict mass-loading factors in star-formation 
driven winds of 0.5-2. We compare these results to those 
of other high-z galaxies and local starbursts in Section 
4.3. 

4.1.3. The Source of the High Velocity Dispersions 



We also investigate the pressure balance in the star- 
forming clumps to determine whether star-formation 
feedback is a viable source of the large velocity disper- 
sions seen throughout the galaxy and by extension, in 
other high-z galaxies. We calculate the contributions 
to the gas pressure from the weight of the gas at the 
midplane (P weight), from star- formation generated tur- 
bulence (Pturb), from radiation pressure (P ra d) and from 
the observed veloc i ty disp ersion (Pkin) from the model of 
lOstriker fc Shettvl (|2011h . hereafter OS11. The model of 
OS11 is based on the assumption that in equilibrium, the 
weight of the gas balances the pressure from feedback. In 
the model, there are two different regimes defined by the 
surface density, where the feedback is either dominated 
by turbulence or radiation pressure; turbulence domi- 
nates in the galactic center from £ flas ~ 100-10 4 M Q pc -2 , 
and radiation pressure takes over for E flas >10 4 M Q pc -2 . 
The pressures are defined as follows: 



Pkin = PO^z 
Pweight = 7rGE 2 /2 

p fp P* v 

4 m* 

Prad = CCKirEEsfr/ 4 



(2) 
(3) 

(4) 

(5) 

where po is the midplane density (which we take as the 
electron density), a z is the observed (narrow-line) ve- 
locity dispersion, E and T,$fr are the surface densities 
of gas and the SFR, f p is a factor of order 1-2 which de- 
scribes the strength of turbulent dissipation, p„/m* is the 
mean momentum injected by SN per unit mass of stars 
formed (which we take as 3000 km s _1 , following OS11), 
e is the mass-to-radiation conversion efficiency from stars 
(which we take as 6.2xl0~ 4 , following OS11) and kjh is 
the mean opacity (which we take as 10 cm 2 g _1 , following 
OS11). 

From these calculations, we find that for clump A, 
Pkin balances P W eight, as proposed by OS11, and that 
Pturfc ro Prad (see Table 3 for details). This latter find- 
ing is consistent with the fact that clump A is close to 
the boundary between the turbulent and radiation pres- 
sure dominated regimes defined by OS11, with E ~ 8000 
M Q pc~ 2 . However, our results are inconsistent with the 
OS11 model in that Pturb+ P rad<P kin (by about an or- 
der of magnitude, even with the errors) , meaning that the 
pressure generated from star-formation feedback, both 
from SN and radiation pressure, does not fully account 
for the observed turbulence in the disk. This result also 
follows from applying equation 2 2 of OS11 to the st ar- 
formation efficiency derived from IGenzel et all ([2010D of 
0.017, wherein the OS11 model predicts an observed ve- 
locity dispersion of 19-37 km s _1 (where the range re- 
flects the the choice of f p ), while we observe dispersions 
of 50-100 km s -1 . 

This issue of the source of the high dispersions 
in high -z galaxies has been expl o red by several au 
thors (IGenzel et all [20(51 l201lt IGreen et al] 
Law et all I2009t lEpinat et alj 120091: IKassin et"atl 



2010: 



2007 



Burkert et al.l l201Q) with no consensus yet reached con- 
cerning the dominant mechanism. Our conclusions on 
this matter are dependent on the exact choice of param- 
eters used, on the uncertainty in our observations, and 
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on the method app lied to derive the intrinsic dispersion 
(jDavies et al.l l2011). However, we are fairly confident in 
our calculation of Pkim as the midplane density that we 
derive from the [SII] doublet is roughly consistent with 
what is observed in dense sta r-forming regions in other 
galaxies (jKewlev et al.ll2001bl ). and several authors have 
found that the high velocity dispersions observed from 
Ha in z ~ 1-2 galaxies is consistent with the disper- 
sions in molecular ga s as observed from CO; for example 
iTacconi et al.l (l2010f); Tacconi et al. in prep for norma l 
SF-galaxies and lTacconi et al.l (|2008t l ; lEngel et al.l (|2010f ) 
for SMGs. 

For clump B, we find slightly different results. As 
for clump A, Pturb+ Prad<Pkm, but unlike for clump 
A, Pkin>Pweight and Pturb>P ra d- We note that the 
first two relations are less significant than for clump 
A since the errors are much higher as a result of the 
more uncertain electron density measurement. The fact 
that Pturb>Prad can be explained by the lower surface 
density of clump B (1200 M Q pc~ 2 ) which falls squarely 
in the turbulence-dominated regime (100-10 4 M pc~ 2 ). 
The results from these calculations support the claim of 
OS11 that the weight of the gas balances the turbulence 
in the disk, which suggests a self-regulatory mechanism 
for SF feedback. However, we find that while SF feed- 
back does contribute to the observed turbulence it cannot 
fully account for it. 

We can look at this question in a different way, by com- 
paring the energy injected into the wind through shocks 
with the turbulent energy in the clumps. To this end, we 
estimate the energy from shocks: 

Eshock — E w f shock, (6) 

where Em is given above and i s hock is the fraction of 
this radiated energy that comes from shocks. Based on 
our analysis in Section 3.3, we choose i s hock values of 0.2 
and 0.3 for clumps A and B, respectively. The turbulent 
energy in the disk is: 

3 

Eturb = 7i M g*> a * ( ? ) 

where M gas is the clump gas mass. We find for the 
E w estimate (E TO = 1/2 M„x v 2 ut ) that ~E shock /E tur b 
is 0.011±°;J 2 4 for clump A and 0.11±g;|| for clump B. 
This suggests that the energy from shocks cannot fully 
account for the energy in turbulence, but in the case of 
clump B, it can contribute significantly. 

This conclusion is consistent with that above, wherein 
SF feedback is an important component to the large 
velocity dispersions, but another strong component is 
necessary. Even the sum of the contributions of star- 
formation feedback (as given by OS11) and shocks falls 
short of the observed energy from turbulence. However, 
if the density of the wind was overestimated as suggested 
in the next section, then the corrected mass and en- 
ergy in the outflow would be ~ 25 times larger, bringing 
^ shock /^turb ~ 1 (as represented by the upper error bar 
for clump B). In addition, the E s hock estimate only in- 
cludes the energy from the warm ionized component in 
the wind, and the hot, neutral and molecular compo- 
nents could also contribute to E shock- Therefore, future 
observations that constrain the contribution to the wind 



from these other components is critical to understanding 
energy balance. 

Assuming our density estimates and mass outflow rates 
are correct, the combination of feedback and shocks can- 
not fully explain the high velocity dispersions observed in 
high-z galaxies. Roughly 85-90% of the velocity disper- 
sion remains unaccounted for and must have some other 
source, such as clump-clump interaction, clump-stream 
interaction, or the release of gravitational energy as gas 
or clumps fall inward. 

4.1.4. The Outflow Can Escape the Disk of the Galaxy 

To test whether or not the outflow can escape the plane 
of the galaxy, we compare the pressures derived from the 
OS11 model with the ram pressure in the wind: 

Pram = PwV out (8) 

where p w is the density in t he wind and is somewhat un- 
certain. iGenzel et al.l ([201 ID adopted a value of 100 cm~ 3 
based on [SII] doublet obser vations of local outflows from 
IHeckman etaLl (|1990l) and iLehnert fc Heckmanl (|1996ft . 
However, more recent work has shown that the actual 
density of the warm ionized component in superwinds 
could be almost 2 orders of magnitude lower (jCecil et al.l 
120011) . Using the larger density estimate (100 cm~ 3 ) for 
clump A, P r am ~ P 'weight, which suggests that the wind 
would barely be able to break out of the plane of the disk, 
while for the lower density estimate (4 cm -3 ), P ram << 
P 'weight- The situation is reversed for clump B which has 
a more evolved and th erefore faster wind (see models by 
IHeckman et al.l (|1990f) which suggest that winds become 
faster with increasing radius) , such that both density es- 
timates yield P ram >> P 'weight- In this case, the wind 
would have no trouble escaping the disk. 

Alternatively, we can estimate the wind density by as- 
suming that Pram balances P weight- With the pressures 
calculated as before, we get n W md — 67 ± 74 cm -3 for 
clump A and 0.4 ± 0.5 cm~ 3 for clump B. This result 
is consistent with clump B having a more extended and 
therefore more diffuse outflow. 

We can also test whether or not gas can escape the 
galaxy by comparing the outflow velocities of the clump 
winds with the escape velocities of the clumps and the 
galaxy. At the distance of the clump HWHM and esti- 
mating the clump mass as 2 x M oas (assuming half of 
the clump mass is in stars), we find that for clumps A 
and B v out /v esc = 0.75 ± 0.39 and 2.4 ± 1.3, respec- 
tively. Therefore, some of the warm ionized gas could 
escape the gravitational potential of the clumps. In fact, 
some of this material could escape the galaxy as well since 
Vout/v e sc,gaiaxy = 0.54 ± 0.26 and 0.99 ± 0.48 for clumps 
A and B. Thus analysis of both the ram pressure of the 
wind and v otlt / v escape suggests that a significant fraction 
of the gas could escape the gravitational potential of the 
clumps distributing the warm ionized gas elsewhere in 
the galaxy and even into the IGM. 

4.1.5. How Clumpy are the Clumps? 

We can estimate the 'dumpiness' of the clumps (or a 
volume filling factor of high density gas) using the ob- 
served electron density (from the [SII] doublet) and the 
calculated surface density. This number can be compared 
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to that for Milky Way GMCs and central starburst re- 
gions to better understand how the star-forming regions 
are distributed within the clumps. We define it as fol- 
lows: 



dumpiness 



b avg 



1 



4/3r c ifim h n H2 



n H2 



rimr 



ii 



T, 



H-2 



(9) 



(10) 



where the latter equation is simple thermal pressure bal- 
ance to obtain the density in the colder molecular gas 
from the hotter ionized gas, and tihii is the local electron 
density derived from the [SII] doublet. We note that this 
thermal pressure balance gives us an upper limit for the 
molecular density since it neglects any non-thermal pres- 
sure due to turbulent cloudlet motion, which would likely 
be more similar for the ionized and molecular phases. A 
volume filling (dumpiness) factor of 1 would mean that 
the high density star-forming gas occupies all of the vol- 
ume of the clumps, whereas a factor of would mean 
it takes up none of the volume. We set T HII = 10 4 K 
and T H . 2 = 30 K (|Elbaz et alJl20ll IHwang et alllMtl 
to obtain ny 2 = 6xl0 5 and 1x10 s cm -3 for the clumps. 
For clumps A and B, we get volume filling factors of 
3.7xl0~ 4 and 2.3xl0~ 4 with la ranges of 0-7.6xl0~ 4 
and 0-6.3xl0~ 4 , respective ly. If, however, we ado pt the 
molecular density found by lDanielson et all (|2011[ ) for a 
lensed z ~ 2.3 SMG of 10 4 cm" 3 , we find much higher 
volume filling factors of 0.022 and 0.0023 for clumps A 
and B, respectively. These different values serve as our 
lower and upper limits for the volume filling factors of 
the clumps and are reflected in the errors (Table 3). 

The volume filling factors based on the local electron 
density of the clumps are almost 2 orders of magni- 
tude smaller than what we find in Milky Way (MW) 
GMCs. Assuming an avera ge density ~ 100 cm~ 3 and 
a local/core density ~ 10 4 (jstahler L Pallal 12005ft . then 
n-avg /niocai (MW) ~ 10~ 2 . Our clump estimates, how- 
ever, are more similar to the volume filling factor in local 
starburst galaxies, based on mod els of M82. In the cloud- 
cluste r starburst region model of lForster Schreiber et al.l 
(2001), they estimate that clouds of molecular material 
are roughly 0.5 pc in radius surrounded by a layer of neu- 
tral and ionized gas up to 1 pc with a spacing of 2-7 pc 
between them. This gives an average volume filling fac- 
tor of 3.6xl0 -3 for the cold molecular gas. This suggests 
that the star forming regions in the clumps are more 
concentrated than regions in bot h MW GMCs and star - 
burst regions unless we adopt the lDanielson et all (|2011f ) 
molecular gas density, in which case, our values are com- 
parable to those from MW GMCs and local starbursts. 

This should not be surprising as the molecular densi- 
ties in the clumps are over an order of magnitude larger 
than what is found in star-forming cores in MW GMCs 
and the pressure in the disk (as seen in 4.1.3 and Ta- 
ble 3) is several orders of magnitude larger than what is 
found in the MW (P ~ 10 -10 ergs cm -1 ) and is more 
comparable to the pressure in local ULIRGs (P ~ 10~ 7 
ergs cm -1 ) (|Downes &: Solomonll 19981 ). Thus, the clumps 
more closely resemble starburst regions in ULIRGs than 
MW GMCs, both in terms of physical properties and spa- 



tial distribution of dense gas. This result implies that the 
very large SFRs of the clumps are not because a higher 
fraction of the volume is involved in star-formation (in 
fact less volume is forming stars than for MW GMCs), 
but because the regions that are forming stars have much 
higher gas fractions. 

4.2. Do the Clumps Represent an Evolutionary 
Sequence? 

In this section, we explore the evolution of the clumps 
through their winds and whether they could typify differ- 
ent stages in clump development and evolution. We have 
seen throughout this paper that while both clumps A and 
B are actively forming stars and produce powerful winds, 
their individual properties are quite different. Clump B 
contains more enhanced shock-like emission line ratios 
and broad emission spatially offset (and therefore fur- 
ther away) from the clump as well as larger line widths, 
both indicative of a more mature outflow. The clumps 
also differ in their volume and surface densities (clump 
A is denser). This could be the result of the expulsion 
of high density gas from star formation feedback from 
clump B or it could be due to the fact that the luminos- 
ity and density measurements from clump B include not 
only clump gas but a large amount of the more diffuse 
gas from the wind. Both explanations further support 
the presence in clump B of a more mature outflow, sug- 
gesting that clump B has been forming stars for a longer 
period of time and thus clump B is older than clump A. 

This notion is supported by the difference in gas 
phase metallicity between the clumps as calculated from 
[Nil] /Ha which suggests an age difference. We calcu- 
late how long it would take clump A to be as metal 
enriched as clump B (assuming the [Nil] /Ha ratio re- 
flects the abundance only and is unaffected by shocks as 
a simplific ation ) .Using th e closed and leaky b ox models 
outlined in|ErB (|2008D (see lGenzel et al.ll201lL Appendix 
C), we find this timescale ranges from 170 Myr - 1.7 
Gyr. If the value of [Nil] /Ha is more strongly affected 
by shocks for clump B (as is suggested by Figure 7) , then 
the resulting timescale would be reduced, but it is likely 
that at least some of the [Nil] /Ha enhancement in clump 
B is due to a higher metallicity, and thus more mature 
clump. 

Given that clump B is likely older than clump A, we 
estimate how long it would take a clump A like object 
to turn into clump B, which has half of the clump A 
M ffas , assuming that the clumps start out with roughly 
the same M sos . Although the gas masses are uncertain 
by 70-85% due to the fact that the physical conditions in 
these ~kpc size clumps could be very different from the 
galaxies from which we derived the KS relation, it is hard 
to imagine that the relative masses between the clumps 
would be that far off, given that the physics and CO to 
H conversion factor in both the clumps are likely similar. 
If all of the extra gas from clump A is turned into stars, 
this would take ~ 200-700 Myr (with the range reflecting 
the SFRs of clumps A and B). However, if outflows are 
removing gas from the clumps as well, then this timescale 
is much shorter; with an r\ of 1, the timescale is twice as 
fast and with an 77 of 3 (corresponding to clump A), the 
timescale is four times as fast. In Section 4.1.4, we found 
that a significant fraction of the outflowing gas could be 
lost from the clumps and perhaps even the galaxy. Based 
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on our finding that rj > 1, the age difference between the 
clumps obtained from the SFR and M gas is at most 350 
Myr, and for this to be consistent with our results based 
on [Nil] /Ha, shocks must play an important role in the 
emission. 

We can also compare clumps A and B with C. Clump C 
has a similar flux and SFR as clump B, and the extended 
[Nil] emission from clump B continues to the region 
around clump C, indicating that it is somewhat enriched 
from star-formation. However, clump C has a narrow 
Ha line profile compared to clumps A and B and has 
no indication of outflows. Perhaps clump C could be 
placed further down the evolutionary sequence of clumps 
A and B: it is declining from its major episode of star 
formation and its outflows are too diffuse to be observed 
or are projected onto clump B due to the inclination of 
the galaxy. 

4.3. Comparison of Star-forming and Outflow 
Properties of ZC'406690 to those of Galaxies at 
Low- and High-z 

High-z SFGs, like ZC406690, in many ways resemble 
both local starburst galaxies and local normal disk galax- 
ies. The results from this work and others suggest that 
the overall structure and star-forming history (steady as 
opposed to bursty) of t hese galaxies more closely resem- 
bles n ormal local SFGs (jWuvts et al.ll2011btlElbaz et al.l 
120111 ) . while properties of the star- forming regions and 
their outflows are more comparable to those of local nu- 
clear starbursts. 

The overall structure and large-scale properties of high- 
z rotation-dominated SFGs are in many ways similar 
to local disk galaxies. First, they both have disk-like 
morphologies (although local galaxies have much thin- 
ner disks) and most of the star-formation in these high-z 
galaxies occurs in an extended disk or ring and not in 
a central and compact starburst region. Second, they 
are not considered extreme star-formers for their red- 
shift (unlike SMGs), and many of these observed galax- 
ies have ki nematic structures sugge s tive of no major in- 
terac t ions dGenzel et aLll2006l. 120081. l201lHShapiro et al l 
20081 iForster Schreiber et al.1 120091: iCresci et al.l 120091: 
Gnerucci et all l2010t Uones et al.l I2010D . Third, with 
the steady accretion of cold gas (in agreement with the 
ACDM model), these high-z galaxies can sustain star 
formation for long timescales (~ 1-2 Gyr). 

In local properties, however, high-z SFGs more closely 
resemble local s tarbursts. From the B PT diagram (Fig - 
ure 6, as w ell as lShaplev et al.l (|2005D : lErb et all (|2006l i: 
ILiu et all(|200ll l we see that normal SFGs at z^l-3 are 
offset from normal local SFGs and better coincide with 
th e region occupied by l ocal starburst galaxies. A s shown 
bv lLiu et all (|2008fl and lBrinchmann et al.l 1)20081 ) . this is 
due in part to variations in ionization parameter, SFR, 
sSFR, and density, between normal SFGs at high and 
low-z, and that in terms of these parameters local star- 
bursts are more similar to t he hig h-z SFGs. It should be 
noted that in the ILiu et al.l (|2008l ) and lBrinchmann et all 
(2008) studies, they compared SDSS galaxies that were 
part of the main star-forming branch and ones that were 
offset from this branch, using similar stellar mass ranges 
for both samples. If we compare the properties (SFR, 
S gas , P, n# 2 , cr, U, volume filling factor) of clumps A 
and B of ZC406690 with local galaxies, we find a closer 



correspondence with starbursts than with normal star- 
forming galaxies (see Table 4). 

The outflows are dictated by the properties of the 
star-forming regions and not the overall structure of the 
galaxy, and are therefore similar to winds observed from 
local starbursts. We find that the mass-loading factor 
for the clumps is consistent with those observed from 
local starburst gala xies: 0.26 for a s ample of several 
IRG s and ULIRGs dRupke et all 12001) and 2.5 for Mrk 
231 (|Rupke fc Veilleuxll2011l)r~Ln addition, fo r Mrk 231, 

Ewind/Lho/ ~ 0.007 (|Rupke fc Veilleuxl l20lTI ). which is 
similar to our estimate for clump A. The measurements 
for Mrk 231 are lower limits; thus if the mass-loading 
and E w j n d/L;, ; factors are larger, then Mrk 231 would be 
more similar to clump B rather than clump A. It is im- 
portant to note that the measurements for Mrk 231 may 
not be directly comparable to the observations presented 
here, as they are made using the NalD absorption fea- 
ture and are therefore sampling a different (and possibly 
more massive and energetic) component of the outflow. 
We also find similar mass loading factors when comparing 
our observations (rj ~ 3-8) to those of wa rm ioniz e d out - 
flows from local dwarf starbursts, where I Martini ([1999D 
finds r] ~ 1-5. 

The outflow properties of individual clumps in 
ZC406690 resemble those from integrated data of other 
high-z SFGs. As for the local starb ursts, most outflows 
from high-z galaxies have rj j> 1 (iPettini et all 120001 : 
iWeiner et al.ll2009tlSteIdel et al.ll2010D . We also find that 
the clumps fit the v m i,t.-M* ajidv ollt -SFR relations found 
at z^l.4 by Wein er et all (|2009| ). In addition, in their 
study of high-z outflows observed in absorpti on of light 
from background galaxies. iSteidel et al.l (|2J110j) find max- 
imum outflow velocities of 700-800 km s _1 , similar to 
what is found for the more rapid clump B outflow (810 
km s _1 ). These similarities suggest that the outflows ob- 
served in previous work from galaxy-integrated spectra 
could originate from massive clumps in those galaxies. 
They also indicate that the star-forming properties in 
high-z galaxies could scale from clump to galaxy size, 
suggesting that the lessons learned from studying indi- 
vidual star-forming clumps could be applied to entire 
galaxies and vice versa. 

If some of the gas can escape the gravitational poten- 
tial of th ese galaxies (as is likely for the hot and fast wind 
fluid, see I Strickland fc Heckm an ( 2009J) and perhaps also 
part of the warm ionized component), these numerous 
and powerful winds arc likely significant in enriching the 
IGM at early times. Such enric hmen t from outflows is 
supported by theoretical models; lErbl (|2008f ) found that 
in their models of metal enrichment by star-formation 
paired with infall and outflow, rj ~ 1 is required to fit 
the mass-metallicity relation at z ~ 2. In addition, the 
results of cosmological simulations by iDave 
indicate that r\ ~ 1 is a necessary condition to reproduce 
the mass-metallicity relatio n as well as IGM abu ndances. 
On the observational side, ISteidel et al.l (|2010l ) find en- 
riched gas as far as 125 kpc away from nearby galaxies 
in a sample of ^100 z^2-3 Lyman-break SFGs. On the 
other hand, if some of the gas falls back onto the clumps 
(as suggested by the fact that v wmd < v esc ,cium P for 
clump A), this could prolong the clump star- formation 
timescales. 
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Evidence for shocks in the clump outflows allows us 
to draw additional comparisons between ZC406690 and 
local and high-z SFGs. The finding that some shock con- 
tribution to the emission is necessary to explain the ob- 
served emission line ratios in the clumps and their winds 
is consistent with observations of outflows from local 
starbursts, which are dominated by shock-like emission 
line ratios even in the absence of AGN. However, if shocks 
are a common feature of high-z outflows, which in turn 
are ubiquitous in z~2-3 galaxies, then this has impor- 
tant consequences f or abundance measur emen ts at high- 
z. As suggested by Shapley et al.l (|2005| ) and lErb et all 
(2006), using [Nil] /Ha alone as an abundance indica- 
tor could be flawed in the presence of shocks or AGN. As 
they note, many of these emission line ratio abundance 
indicators are calibrated from galaxies and star-forming 
regions that populate the HII branch of the BPT dia- 
gram. Therefore, galaxies that are offset from this branch 
(mainly in terms of [Nil] /Ha ), such as local starbursts 
and high-z SFGs, require a different calibration. An- 
other option is to use photoionization and shock models 
with multiple emission line diagnostics to get more accu- 
rate abundances. It is not yet clear whether abundances 
are overestimated due to shocks in most observed high- 
z SFGs, but this is an issue that deems more attention 
through continued high-resolution AO emission line mea- 
surements in these systems. 

5. CONCLUSIONS 

Using high S/N SINFONI/VLT observations with AO, 
we find powerful outflows localized to the star-forming 
~kpc-size clumps of the z ~ 2 SFG, ZC406690. Spectra 
of the clumps and surrounding regions have broad, blue 
Ha line wings extending up to 800 km s _1 from line cen- 
ter. Much of this broad emission is coming from a region 
offset from the center of clump B that also has enhanced 
[Nil] /Ha and [Sll]/Ha ratios, indicative of a superwind 
laden with shocks. The outflow velocities derived from 
the line widths are comparable to and slightly in excess 
of the escape velocity of the clumps (and the galaxy), 
indicating that some of the wind material could fall back 
on other parts of the galaxy or perhaps escape out to the 
IGM. 

By comparing the en ergy and momentum injection 
rates of the wind to the iMurrav et al.l (|2005| ) model of 
a starburst-driven wind, we find that the clump A wind 
satisfies the predictions for both energy- and momentum- 
driven winds, but that the clump B wind cannot be ex- 
plained by either within the context of this model. Star 
formation feedback from the winds could greatly con- 
tribute to the large velocity dispersions seen in ZC406690 
and other z ~ 2-3 galaxies, although based on measure- 
ments from the clumps, feedback alone cannot account 
for this turbulence. 



The clumps have a large range in ages (from metallic- 
ity), which is reflected in the maturity of their outflows 
(as seen in the variation in outflow velocity, physical ex- 
tent, mass- loading and shock- like emission line ratios), 
and they appear to be more extincted than the surround- 
ing regions, as one would expect for dense star-forming 
regions. It is possible that the younger clump (A) re- 
sembles what clump B looked like at an earlier time and 
that the current state of the galaxy gives us a snapshot 
of various stages of clump formation and evolution. 

In terms of global properties (disk-like morphology, 
location of star- formation, star- formation timescale), 
ZC406690 is similar to local disk galaxies, while in terms 
of local properties of the star-forming regions and the 
outflows (SFR, SsFfi, nf/ 2 , pressure, ionization parame- 
ter, Vout, mass- loading, volume filling factor), its clumps 
more closely resemble star-forming regions in local star- 
burst galaxies. The outflow properties of the clumps are 
also similar to those of outflows observed in other galax- 
ies at z ~ 1-3, suggesting that massive clumps could be 
the source of a significant fraction of outflowing material 
in these galaxies as well, and that although ZC406690 is 
quite massive, it may not be very unusual for z ~ 2 SFGs 
in terms of star-formation feedback. 

These types of winds coming from kpc-size regions 
could have a significant effect on galaxy evolution dur- 
ing this epoch. They could inhibit star-formation with 
the injection of heat into the ISM and the mass-loss as- 
sociated with the outflowing material in the wind could 
reduce the lifetime of the clumps. They also influence 
the chemical composition of the ISM and the IGM by 
expelling metal-rich gas from star-forming regions. The 
results presented in this paper inform our understanding 
of the properties, evolution, and feedback of star- forming 
regions at high-z, which in turn teaches us about the role 
that these star-forming regions play in galaxy evolution 
on Gyr timescales. 
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TABLE 1 

Emission Line Data from the Clumps 



Region 


% broad 


C" narrow 


®broad 


[Nil] /Ha 


[Sll]/Ha 


% shock" 


x^fit 


Clump A 


50 ± 4.8 


82.7 ± 4.2 


199 ±11.9 


0.11 ± 0.014 


0.14 ± 0.029 


5 ± 5 


0.78 


Wind A 


44 ± 7.2 


84.7 ± 4.9 


197 ± 20.5 


0.15 ± 0.036 


0.21 ± 0.061 


10 


0.64 


Clump B 


46 ± 4.6 


81.9 ± 4.6 


244 ± 68.2 


0.20 ± 0.031 


0.20 ± 0.055 


15 + ' & 
10 -15 


0.59 


Wind B 


71 ± 2.9 


97.7 ± 6.0 


291 ± 12.7 


0.32 ± 0.032 


0.27 ± 0.068 


30 +3U 
JU -20 


0.58 



a The shock contribution is calculated by comparing the data to Figure 10 of lRich et al.l (2011). 
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TABLE 2 

Observed Properties of the Clumps 





Clump A 


Clump B 


E(B-V) a (Ha/H/3) 


0.493 ± 0.165 


0.493± 0.165 


E(B-V) (SED) 


0.33 ± 0.1 


0.33 ± 0.1 


[Oll]A3729/3726 


1.38 ± 0.41 


1.51 ± 1.14 


n e ([Oil]) (cm- 3 )" 


80, [l-100] c 


5, [1-1000] c 


[SlI]A6716/6731 d 


0.7 ± 0.1 


1.2 ± 0.3 


n e ([SII]) K 3 )W 


1800, [1100-3000] c 


290, [30-750] c 


[OlII]/H8 a 


6.72 ± 1.31 


4.62 ± 10.30 


[Nil] /Ha c 


0.11 ± 0.018 


0.21 ± 0.029 


[SII] /Ha e 


0.18 ± 0.049 


0.18 ± 0.067 


[OlI]/[OHI] e 


0.99 ± 0.14 


1.06 ± 0.35 


U (from O32) (cm- 1 )"^ 


5.50xl0 7 , [4.75xl0 7 -l.lxl0 8 ] c 


8.75xl0 7 , [3.75xl0 7 -1.2xl0 8 ] c 


U (from O32) (cm- 1 )''' 




1.03x10 s , [4.50xl0 7 -1.25xl0 8 ] c 


SFR (M s yr- 1 )' 


40 ± 12 


11 ± 3.3 


r ut,HWHM (kpc)-f 


0.8 ± 0.3 


1.2 ± 0.5 


U (from SFR) (cms" 1 ) 


7.4xl0 7 ± 6.1xl0 7 


8.3xl0 7 ± 1.0x10" 


12 + logl0(O/H)' 1 


8.42 ± 0.12 


8.62 ± 0.096 



a For clump B, these numbers are estimated assuming Ha/H/3 is the same for clumps A and B. 
b n e is calculated from Ostcrbrock (1989) assuming collisional de-excitation with T e = 10 4 K. 
c la range 

d from lGenzel et al.l S20TH) 

e For regions including both the clump and the wind which are spatially consistent among H, J, and K-band cubes. 

f For clump B, assuming Ha/H/3 is the same as for the diffuse gas. 

9 For this calculation, metallicities are assumed to be those as calculated using the combined method of Kcwlcy & Dopita (2002). 
h Using the [Nil] /Ha conversion of lDenicolo et al.l (f200W) . 



APPENDIX 

APPENDIX A: SYSTEMATIC ERRORS FOR M GAS 

In this section, we derive the systematic errors in our M 9as calculation from a combination of (1) scatter in the 
intrinsic relation, (2) the varying slope and zero point of the relation as quoted by different authors, (3) the extrapolation 
of the relation from low-z to high-z, and (4) the application of the relation on ^kpc scales. 

First, we consider the scatter of T, gas within an individual Kennicutt-Schmidt (KS) relation. iBigiel et al.l ([201 1) 
compare T, gas and T,sfr for thousa nds of positions in 3 nearby disk galaxies, and find that the la scatter is roughly 
0.3 dex or 70% (see their Figure 1). Ge nzel et al.l (|2010D find a similar systematic uncertainty for T, gas of 0.3 dex in 
the KS relation for both local and high-z galaxies. 

Second, we consider th e syste mati c uncertainty r esulti ng from extrapolating the low-z KS relation to a high-z galaxy. 
As seen in iGenzel et al.l (|2010f) and Daddi et al. (2010), high and low-z galaxies on the star- forming main sequence 
(that is, neglecting local ULIRGs and high-z SMGs) fit on the same KS relation with an offset in zero point of less 
than a few x 0.1 dex, in the sense that high-z galaxies could have slightly lower gas masses than the local relation 
would predict. ZC406690 falls on the upper SFR and M* end of the SF main-sequence. 

Third , we consider the scatter between different versions of the KS r elation by looking at the most recent observational 
res ults (IGenzel et al.ll2010l ; iDaddi et al.l l2010; Kc nnicutt et al.1 12007). We do not consider the slope and zero-point of 
the IBigiel et al.l (|201l( T relation, as they note that theirs is not a rigor ous fit. From the thre e other r elations, the 
maxim um difference in T, gas for our clumps is around 50%. Both the iKennicutt et al.l (|2007l ) and the IDaddi et al.l 
(2010) relations have slightly larger slopes, so the effect of using these would be to lower the resulting gas masses. 
We note that the gas masses rely on the extinction correction based on the galaxy-integrated SED fitting and not 
on the Ha/H/3 measurements (which are uncertain for clump B). As mentioned in section 3.1, applying the Balmcr 
decrement-based extinction correction would increase the gas masses by 50%, the opposite effect of using the other KS 
laws. 

The statements made above rely on the KS conversion for normal main sequence star-forming galaxies. While 
ZC406690 is in fact on the star-forming main sequence, there is some indication that the clumps might follow the 
ULIRG KS law, both because the clumps have similar properties to local starburst regions and also based on recent 
work by ISargent et al.l (|2012f ) suggesting that z ~ 2 normal SFGs can simultaneously have a normal star-forming 
component and a starburst component. Using t he ULIRG law to determine the gas masses for our clumps would lower 
the resulting masses. From IGenzel et all (|2010l) , Fi gure 4, we find that the offset from the main-sequence KS law to 
the ULIRG law is around 0.8 dex for Hsfrs relevant for the clumps. Therefore, we assume a lower limit for the gas 
masses of 0.8 dex. 
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TABLE 3 

Calculated Properties of the Clumps and their Outflows 





Clump A 


Clump B 


SFR (M jo— 1 )" 


40 ± 12 


11 ± 3.3 


M ou t(lvl0 yr ± ) Ll 


-1 1 7 +111(0 

117 -95 


7fl +b90 
'8 _ 78 


v out (km s _1 ) a 


440 ±150 


810 ± 270 


Vout /"V escape, clump 


0.75 ± 0.39 


2.4 ± 1.3 


^ oui f~^e scape , a alaxy 


0.54 ± 0.26 


0.99 ± 0.48 


trf,™ (Myr) 


1.8 ± 0.9 


1.5 ± 0.8 


M U ,(M ) 


1 79x10 s + 1£,6 >< loa 


1 09x10 s + 8 - 36 xi0° 

ij . 4t :} A 111 


Eu,(ergs) 


3 46x10 56 +3-81x10°' 

0.40X1U -3.21X10 56 


7 IfixlO 56 +5-50X10 5 ' 


m , rivr,~+ 

aL gas,mol \^LQ ) 


l.bxlU „ 13xl0 io 


7 8xl0 9 + 5 ' 5 xlO a 
rsxiu _ 66xlo9 


S 9 a S (M Q pC- 2 ) a 


„„„„ +6000 

8000 telZ 


1 nnn +840 

1200 t m > Q 


Ssfjj (Mopc-^yr- 1 ) 11 


13.5 ± 4.0 


1.6 ± 0.5 


Pfein(dyne cm" 2 ) 


2.9xl0" 7 ± 1.6xl0~ 7 


5.3xl0~ 8 ± 5.5xl0" 8 


Plight (dyne cm -2 ) 


3 0x10~ 7 + 2 - 9x10 
d.UxlU _ 36xl0 -7 


6 7x10 -9 + 6 ' 6x l0 


Pturb(dync cm" 2 ) 


6.8xl0~ 9 (± 2.1xl0~ 9 ) x f p 


8.4x10 (± 2.5x10 ) x i p 


Prod (dyne cm ) 


'• ixiU -6.4x10-9 


1 3x10-1° +1.0X10- 1 " 
1.3x10 2xio-io 


Pram (dyne cm -2 ) 


4.4xl0- 7 (+^^°_ 7 )x( Pu ,/100) 


l.sxio- 6 (if°*;°_ 6 )xowioo) 


Punnd (cm -3 )' 


67.3 ± 73.7 


0.4 ± 0.5 


Lfcoi (ergss -1 ) 


1.6xl0 4b ± 4.8xl0 44 


4.4xl0 44 ±1.3xl0 44 


M out X Vout 2 / L boi 


0.0045, [7.64xl0" 5 -0.048] c 


0.046, [0.00021-0.33] c 


Pout/Prad 


6.2, [0.5-65. l] c 


33.7, [2.1-244.5] c 


r? (M out /SFR) a 


2.9 


8.6 + b 7 %' 


E s „ oc fe (ergs) 


3 46 x10 55 +3 8lxlO :>D 
6.10 X1U _ 3 . 65xl0 55 


2 1t; x in56 + 1. 66X10 1 " 
^.l&xlU -3.86X10=6 


E tur6 (ergs) 


3 28x10 57 + 2 -31xl0°' 
d.^8xlU _ 2 .77 X1() 57 


1 89xl0 57 +1-33X10 3 ' 

i.syxiu _ 159xl0 57 


£» (n a u 9 /ni oca ;) 


2.2xl0" 2 -3.7xl0" 4 


2.3xl0- 3 -2.2xl0- 4 



a from lGenzel et al.l 120TH ) 

b assuming P ram ~ Plight 
c la range 



TABLE 4 

Comparison of Clump and Galaxy Properties 





Clumps 


ZC406690 


Milky Way 


ULIRGs 


SFR (M yr" 1 ) 


10-40 


300 


3 


100-1000 


S 9as (M pc- 2 ) 


1200-8000 


1100 


200 


3000-30,000 


f 9 as 


0.4-0.7 


0.3 


0.1 


0.1 


logl0(P/k) a 


8.6-8.9 


9.0 


5.9 


9.3 


logl0(nH,,) 


4.0-5.8 




2 


3-4 


a (km s — 1 ) 


85 


63 


5 


50-100 


loglO(U) (cms" 1 ) 


7.8-8 




7.5 


7.7 




3-8 






0.1-2.5 


P/(Uol/c) b 


6-33 






0.1-0.5 


dti /t b 

— / L iR 


0.04-0.005 






0.0001-0.007 




2xl0~ 2 -4xl0- 4 




10~ 2 


3.6xl0~ 3 



a from Pfcin 

'These values are derived from Ha for ZC406690 and from Na ID for ULIRGs URupke et al.ll200"5t IRupke fc Ve illcux 2012). However, this 
might not be a fair comparison since these measurements sample very different components of the wind. In addition, the largest values for 
ULIRGs (for -q and ^f/L/ fl ), which come from Mrk231, are lower limits ijRupke fc Veilleux|[20Tll ) . 



Fou rth, we believe it is reasonable to assume the global KS relation for regions of ~1 kpc in size. iBolatto et al.l 
(2008), fi nd that many of t he gas properties of extra-galactic GMCs are consistent with properties of the MW. In 
addition, iBigiel et al.l ( 20111) find that their b inned 1 kpc-scale relation is almost identical to the relation binned from 
global literature data. iKennicutt et al.l (|2007l ) also look at th e KS relation on t he scale of 0.5 and 2 kpc in M51a and 
find a very similar slope to that from the global relation from iKennicuttl (|1998[ ). and an offset in zero-point of around 
0.4 dex. It should be noted that the M51a integrated point falls 0.24 dex below the global law, so a good estimate for 
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the offset in zero-point from extrapolation of the global to local KS relation would be a roun d 0.15 dex. In addition , 
recent observations of the lensed cosmic eyelash galaxy at z~2.3 bv lSwinbank et all (|201lD and Danielso n et al.l (|2011h . 
probe the KS relation on the scale of 100 pc submillimeter-bright clumps. They find that the global KS relation should 
hold on scales of roughly ^kpc, but should break down below 100 pc, and that globally, the galaxy is consistent with 
local ULIRGs (as expected for a high-z SMG). Finally. iLerov et all (|2012D study Ha-, UV-, and IR-derived SFRs in 
addition to IRAM CO data at 1 kpc resolution for 30 disk galaxies, and find that the Ha-derived "Esfrs could be off 
by about 0.3 dex at small vs. large scales, comparable to the 0.3 dex systematic uncertainty we assume. 

The systematic uncertainties from applying the local KS relation to high-z galaxies, from the scatter between 
relations, from extrapola ting to ~kpc scale s and from propagating t he T<sfr error are all less than or equal to the 
intrinsic scatter from the iBigiel et al.l ([201 If ) and iGenzel et al.l (|2010f ) relations of 0.3 dex, so we therefore use this as 



the upper limit of the systematic uncertainty for our derived clump Y, gas and as mentioned above, we use 0.8 dex as 
a lower limit. 

APPENDIX B: DERIVATION OF MASS OUTFLOW RATES 

In this section we discuss the m o del u sed to derive the mass outflow rates of the clumps, which is identical to 
the derivation used in IGenzel et all ()201lD . We derive the mass in the outflow and the mass outflow rates assuming 
two different geometric models. For both, we assume an outflow of warm, ionized gas into a solid angle SI with 
radially constant outflow rate and velocity. As we are only calculating the mass outflow rate for the warm-ionized 
component and neglecting other components of the wind (neutral, molecular, and hot), our masses and rates are 
therefore lower limits. We assume case B recombination of photoionized gas with T = 10 4 K, giving a volume emissivity 
of 7fl Q = 3.56 x 10~ 25 erg cm~ 3 s _1 . We then calculate the mass of outflowing material from the extinction-corrected 
Ha luminosity: 

L Ha = lHa J nR 2 n e {R)n p (R) dR, (Bl) 
M HII He = 1.36 x m p [ nR 2 n p dR = lM *rn p L Ha ^ (B2) 

J jHan e ff 

M out = fLR 2 1.36 x m p n(R)v out = M H n,He X ^ (B3) 

where L# a is the extinction-corrected Ha luminosity, v out is the outflow velocity derived from the Ha line profiles 
and given in Table 3, and Rout is the maximum radial extent of the outflow, which is model dependent, as is n e // 
(electron density) . The extinction-corrected Ha luminosity is taken from the broad-component of the line only, and is 
multiplied by an additional factor of 2 to account for the asymmetry of the broad line wing (i.e. much of the red wing 
is obscured, see Figure 2). 

For our first model, we assume that the average electron density in the outflow scales like R -2 , but the local electron 
density of clouds and dense filaments in the outflow (from which we observe Ha emission) is constant at n e = 100 
cm~ 3 . In this model, R ou t ~ Rhwhm and is given in Table 2. As mentioned in section 4.1.4, the density in the wind 
could be closer to 4 cm -3 and this is reflected in the uncertainties for M out . In fact, the value of n e // is the dominant 

source of uncertainty for M„, and M out and for subsequent calculations that depend on these quantities. We find that 
the lower limit for n e ff (10 cm -3 ) gives unrealistic upper limits for the mass-loading factor of 30-60 (see Table 3). 

For the second model, we assume that the entire outflow cone is filled with the ionized, Ha-emitting gas, and that 
n e // still scales like R~ 2 . Thus, n-eff — rii n Ri n j R ou t, where Ri„ is the radius at which the out flow is launched an d 
n in is the corresponding electron density at that radius. Assuming the outflow is very extended dSteidel et al.l feOlO'). 
we take R ou t ~ Rdisk ~ 10i?j n . From the n 2 scaling of the Ha emission, Rhwhm ~ 2.3iij n . For nj„ we take an 
approximate average of the clump electron densities derived from the [SII] doublet, ~ 1000 cm -3 . Thus, for the 
second model, R ou t = 4.3 kpc and n. e ff = 100 cm~ 3 . These two models likely span the lower and upper limits for R OM t 
and so we take the average of the results from the two models as our wind mass and mass outflow rates (see Table 3). 



